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I. 



We present numerical results for the magnetic susceptibility of a Heisenberg 
antiferromagnetic spin ladder, as a function of temperature and the spin-spin 



£h , interaction strengths J± and Jj|. These are contrasted with new bulk limit 

results for the dimer chain. A fit to the experimental susceptibility of the 

candidate spin-ladder compound vanadyl pyrophosphate, (VO)2P207, gives 

the parameters Jj_ = 7.82 meV and Jii = 7.76 meV. With these values we 

predict a singlet-triplet energy gap of E gap = 3.9 meV, and give a numerical 

estimate of the ladder triplet dispersion relation u(k). In contrast, a fit to 

the dimer chain model leads to J\ = 11.11 meV and Ji = 8.02 meV, which 

predicts a gap of E gap = 4.9 meV. 
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I. Introduction 

Quantum antiferromagnetism in lower- dimensional systems has proven to be a very rich 
subject. Some of the more dramatic developments include the realization that half-integral 
and integral spin chains have very different excitation spectra, PJ and evidence that two- 
dimensional antiferromagnetism is a crucial component of high temperature superconduc- 
tivity. 0-g 

The Heisenberg spin ladder is interesting theoretically as an intermediary between half- 
integer (S = 1/2) and integer (S = 0,1) spin chains. This system has isotropic nearest- 
neighbor interactions along the chains (Jm) and along the rungs (J±) of a ladder geometry, 

H = J||£S i -S i + J±Y,Si-Sj • (1) 

~ I 

Previous studies of the Heisenberg spin ladder have discussed the ground state energy 
and the dependence of the singlet-triplet energy gap on J±/J\\ = on in the antiferromag- 
netic regime, |5],|7j the triplet spin-wave dispersion relation, the behavior of the system 
under doping || and the dynamical structure function. |7j Ferromagnetic rung couplings 
(J_i_ < 0, J|| > 0) have been studied using Lanczos and Monte Carlo techniques and the 
renormalization group; |J these references suggest that a gap exists for all Jj_ < 0. 

In Ref. @ we presented numerical and analytical results for the ground state energy and 
triplet spin- wave dispersion relation of an S = 1/2 spin ladder, as well as numerical results for 
the structure function S(k,ui). We found evidence that a singlet-triplet energy gap appears 
for any interchain coupling J±/ J\\ > 0, and that the spin-wave band minimum is at k = it, 
but the band is folded so the maximum energy occurs between k = tt/2 (for J±/J\\ = 0) and 
k = (for J±/J\\ = oo). The bandwidth was found to be relatively insensitive to the rung 
coupling J_|_, and varied between ttJ\\/2 (for J±/J\\ = 0) and 2J\\ (for J±/J\\ = oo). 

The antiferromagnetic spin ladder may be realized in nature by the antiferromagnet 
vanadyl pyrophosphate, f|,(| (VO) 2 P 2 7 . This material has a ladder configuration of S = 
1/2 V +4 ions (Fig. 2 of Ref. 0), with spacings of 3.19(1) A between rung ions and 3.864(2) 
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A between chain ions, and has a magnetic susceptibility characteristic of an antiferromagnet 
with an energy scale (from the susceptibility maximum) of ~ 7 meV. The closely related 
material V0(HP04) -| H2O has isolated V +4 ion pairs at a similar separation of 3.10 A, 
and is well described magnetically by independent spin- 1/2 Heisenberg pairs with a coupling 
(in our conventions) of J = 7.81 meV. |T0[ 



Although the (VO)2P2C>7 lattice clearly shows a ladder configuration of V +4 ions, they 
might in principle interact magnetically as a different spin system. This was the case for 
Cu(NOs)2-| H 2 (Ref. (TT|]), which was originally considered to be a possible spin ladder 
system but was subsequently found to interact as a dimer chain, described by the Hamilto- 
nian 

H = 7^ < J\ §2i • S2i+1 + J2 ^>2i+l ' ^>2i+2 \ i (2) 

i 

with Jil J\ = Oid (0 < ad < 1). It was not practical to distinguish between ladder and dimer 
chain models of copper nitrate from bulk thermodynamic properties alone, which were found 
to be very similar for the two systems. [|lTJ The issue was finally decided in favor of the dimer 



chain model by proton resonance |12| and neutron diffraction ]T3[ experiments. 



Vanadyl phosphate presents similar ambiguities. Although its susceptibility has been 
measured and is accurately described by the susceptibility of a dimer chain [[| (2), it is 
widely believed that the ladder Hamiltonian (1) will lead to a very similar x{T) and so 
is not excluded by the good agreement with the dimer chain. Since no theoretical results 
have been published for the bulk limit ladder susceptibility in the relevant Jj_ ~ Jh regime, 
comparison of the experimental susceptibility to the ladder model has not been possible. In 
this paper we present new numerical results for the bulk limit susceptibility of ladders and 
dimer chains, and fit these to the data for (VO)2P2C>7 . As we shall see, these two models 
do give very similar results for the susceptibility, and both give excellent fits to (VO)2P2C>7 
with appropriate parameters. The ladder x{T) is preferred, although the differences may be 
less important than the approximations made in the models. 
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II. Spin ladder and dimer chain susceptibilities 

We determined the susceptibility on finite lattices by generating all energy levels {Ei} and 
their multiplicities {d{\ in each sector of fixed total S z , using a Householder algorithm. The 
susceptibility was then obtained through its relation to the expected squared magnetization, 
summed over energy levels and total S z sectors; 

frr\ 2 2 o l^S z Zji ° z a i e / \ 

X(T> = 9 ^ Ei * e-*e. ■ (3) 
This approach has the advantage that explicit eigenvectors are not required. For the ladder 
geometry we used (3) to determine \ on 2 x L lattices with L = 3, 4, . . . , 8 for couplings 
J±/ J\\ — a i — 0-5, 0.7, 0.9, 1.0 and 1.1 and for a range of T / Ju values; for reference purposes 
the 2x8 results are given in Table I. To estimate the bulk limit susceptibility the 2 x L 
results were fitted to the form \l = Xoo + ae~~ bL at each T and coupling, independently 
for L = even and L = odd, since these approached the bulk limit from opposite directions. 
This gave independent estimates of the bulk limit susceptibility, which allowed a test of the 
accuracy of our extrapolation in L. Our bulk limit estimate was taken to be the average 
of the even-L and odd-L extrapolated values. For the dimer chain we followed a similar 
procedure for parameter values Jij J\ = ad = 0.2, 0.4, . . . , 1.0 and L — 4, 6, . . . , 16; Table II 
gives the L = 16 results. For the dimer chain there is no odd/even effect in the long axis (we 
always assume an even number of spins so the ground state has no net magnetization), so we 
had only a single extrapolation in L for each T. This was compensated by smaller finite size 
artifacts than on the ladder, because the long axis of the dimer chain spanned a maximum 
of 16 rather than 8 spins. For fitting purposes it is useful to have a parametrization of these 
results that allows accurate interpolation in T and interaction strengths. We tested several 
forms and found that the six-parameter function 

X (T) = | [1 + (T/e 2 r (e c * /T - 1) ^ [1 + (c 5 /T) Ce ]^ (4) 

adequately describes both the experimental data JEJJ and the theoretical ladder and dimer 
chain susceptibilities over a range of parameters relevant to (VO)2P207 . This form also has 
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the advantage that it incorporates the exponential behavior expected at low temperatures, 
unlike other parametrizations used previously for the dimer chain, and at high temperature 
it gives the correct Curie form g 2 fi 2 B /4kBT. (The overall normalization c\ is identically equal 
to g 2 fij 3 /4:kB in all cases, but the V 4+ ion g- factor is unknown a priori and is determined 
when c\ is fitted to the data.) 

For the two theoretical susceptibilities we fitted our numerical bulk limit results to (4), 
with each of the coefficients c 2 , . . . , cq taken to be quadratic in the ratio of the two coupling 
constants. Preliminary fits to the (VO) 2 P20 7 susceptibility indicated that the values at = 
1.0 and ad = 0.7 were close to optimum, so we parametrized our bulk limit results in terms 
of the departure from these values. The fitted coefficients c 2 . . . c 6 for the ladder over the 
range 0.9 ^ ^ 1.1 (with J\\ = 1) were found to be 



c 2 


= +2.315- 


4.035(c^ - 


1.0) + 7.050(a<- 


1.0) 2 


(5a) 


c 3 


= +0.403 - 


1.025(a< - 


1.0) + 1.850(a<- 


1.0) 2 


(5b) 


c 4 


= +0.443 + 


0.225(a< - 


1.0) + 0.850(a*- 


1.0) 2 


(5c) 


c 5 


= +0.745 - 


0.390(a< - 


1.0) + 0.200(a<- 


1.0) 2 


(5d) 


c 6 


= +1.628 - 


0.110(a< - 


1.0) + 0.600(^- 


1.0) 2 


(5e) 



and for the dimer chain (with J\ = 1) over the range 0.6 ^ a d ^ 0.8 we found 



c 2 


= +8.145 - 61.76(a d 


- 0.7) +406.3(a d 


-0.7) 2 


(6a) 


c 3 


= +0.562 + 2M0(a d 


- 0.7) + 7.300(a d 


-0.7) 2 


(6b) 


c 4 


= +0.456 - 0.435(a d 


- 0.7) + 1.750(a d 


-0.7) 2 


(6c) 


c 5 


= +0.592 + 1.090(a d 


- 0.7) + 0.400(a d 


-0.7) 2 


(6d) 


c 6 


= +1.663 + 0.160(a d - 


-0.7) - 1.40(a d - 


0.7) 2 . 


(6e) 



Due to the presence of large coefficients this parametrization is not useful far from the 
parameter ranges cited; if required the coefficients could be determined directly from the 
bulk-limit numerical results. 
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Our numerical results for the extrapolated bulk limit susceptibility of the ladder and 
dimer chain are shown in Figs.l and 2 respectively, together with the interpolating functions 
defined by (4-6). The interpolating functions reproduce the bulk limit susceptibility with a 
typical accuracy of a few times 10 -4 over the parameter ranges quoted above. 

III. Comparison with the experimental (VO) 2 P20 7 susceptibility 

In a previous study Johnston et al. |§ presented results for the susceptibility of 
(VO) 2 P 2 7 , and noted that the susceptibility of a dimer spin chain gives a very good 
description of the data for ad = Jij J\ = 0.7 (fixed from an interpolation of theoreti- 
cal curves for 0.6 and 0.8 from Refs. J\ = 11.32 meV (hence J 2 = 7.93 meV) and 
g = 2.00. (Note that in our conventions the {J n } are twice as large as in Refs. || and ||.) 
A similar coupling of J = 7.81 meV was determined for VO^PO^- 1 ; H 2 0, which consists 
of isolated V +4 dimers, |K|] and the g factor of the V +4 ion is known to be quite close to 



2 from studies of other vanadium phosphates. |15j Of course it is not clear how the fitted 
dimer chain parameters relate to (VO) 2 P 2 07 if it proves to be a spin ladder. 

To confirm these results we fitted our three-parameter (Ji,ced,g) dimer chain suscepti- 
bility, described by (4) and (6), to the data of Ref. ||, which consists of 606 values of x{T) 
from T = 7.2°K to 344.34°K. We found the optimum parameter values to be J\ = 11.11 
meV, ad = 0.722 and g = 1.99. These are essentially the parameters found by Johnston et 
al., and the minor differences are presumably due to the systematic errors in interpolation 
(perhaps 1% in parameter values). The fitted dimer chain susceptibility and the data for 
(VO) 2 P 2 07 are shown in Fig. 3. 

We similarly fitted the three-parameter (J||, a^g) ladder susceptibility (4), (5) to the 
experimental (VO) 2 P 2 07 x(T) data over the full temperature range. The optimum ladder 
parameters were found to be 

J,, = 7.76 meV , (7a) 
a t = JJJ\\ = 1.007 , (7b) 
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g = 2.03 . (7c) 

The proximity of g to 2 provides a plausibility test of the fit, as does the fitted value of 
J_i_ = 7.82 meV, which is almost identical to the isolated-dimer J = 7.81 meV found previ- 
ously [0 in VO(HP0 4 )-| H 2 . These results suggest that (VO) 2 P20 7 is very close to a 
uniform ladder (J± = Jm), which is presumably accidental because the rungs are bridged by 
two oxygens, whereas the chains have only single oxygens between V +4 ions. The fitted lad- 
der x{T) is shown in Fig.4, and this model evidently also gives a very good description of the 
experimental data. The goodness of fit, defined by the residual / = J2i(Xex P t.(Ti)— Xthy.iTi)) 2 1 
slightly favors the ladder model over the dimer chain. We cannot choose between the models 
definitively from the susceptibility data, however, because the variation in (VO) 2 P20 7 sus- 
ceptibility estimated from samples with different annealing histories (Fig. lb of Ref.( H])) is 
somewhat larger than the difference between the predictions of the ladder and dimer chain 
models. 

Finally, to test how well J± and J\\ are determined, we studied the residual / in con- 
strained two-parameter fits with g and J\\ variable but at = J±/J\\ fixed. As we changed at 
from the optimum value 1.007 we found that by 0.90 and 1.12 the residual had increased by 
a factor of two. As we increase ae through the range [0.9, 1.1] the fitted value of Jm decreases 
from 8.2 meV to 7.3 meV, which can be taken as a conservative estimate of the accuracy to 
which J\\ is determined by the susceptibility. The fitted g factor remains close to 2.03 over 
this range. Outside this range of at there is a rapid decrease in the quality of fit, reaching 
a factor of five increase in / by at = 0.81 and 1.24. 

IV. Predictions of the ladder and dimer models 

Since we have determined ladder parameters for (VO)2P2C>7 from our fit to the sus- 
ceptibility, we can use the results of Ref. |7| to give predictions for the gap and spin-wave 
excitation spectrum. From Fig.2 of that reference we can see that the gap near J±/J\\ = 1 
is quite well determined by the Lanczos and Monte Carlo studies. An approximate linear 
interpolation gives 

7 



E, 




0.50(1) + 0.65 (a t - 1) , 



(8) 



so for the optimum fitted parameters we predict 



E 



gap 



3.9(1) meV . 



(9) 



Over the parameter range 0.9 < J±/ J\\ < 11 discussed above the predicted gap increases 
from 3.6 meV to 4.1 meV. |TjJ 

As was noted in Fig. 3 of Ref. M, for J± = J\\ the minimum energy required to excite a 



triplet spin wave on the ladder as a function of k closely resembles the dispersion relation 
of a spin-1/2 chain, except for the presence of excitation gaps. The lowest excitation is 
at k — 7r, where the gap is uj(tt) = 0.50Jn, the maximum is shifted to a k < 7r/2, and a 
secondary minimum is at k = 0. This dispersion relation is symmetric about k = it. 

A complication not noted in Ref. |7j is that the lowest-lying triplet spin-waves with these 
parameters arise from two distinct bands. The "primary" band, which contains the lowest 
gap, is odd under chain interchange (k± = it), and is shown as solid lines for J±/ Ju = 0.5, 1.0 
and 2.0 in Fig.5. For large J± these are excitations of a single rung, with energy u « J±. 
The "secondary" band (dashed lines in Fig.5) is even under chain interchange, and for large 
J_i_ these states consist of two excited rungs (hence u' ~ 2J± and the even symmetry), 
with the two 5 = 1 excited rungs coupled to give S tot = 1. Thus the secondary band 
may be interpreted as the excitation of two spin-wave quanta of the primary band. This 
interpretation leads us to anticipate several features of the secondary dispersion relation in 
the bulk limit, for example u'(k = 0) = 2u)(k = tt), so the band minimum of the secondary 
band in (VO)2P2C>7 should lie at 7.8(2) meV given our parameters. One may similarly 
construct the entire secondary uj'(k') given the primary u(k) (assuming there are no bound 
states), by finding the minimum-energy combination of two quanta with specified k! . 

In our representation in Fig.5 we fitted the function 



(k) = [u{0) 2 cos 2 (A;/2) + ^(tt) 2 sin 2 (fc/2) + c 2 sin(fc) 2 ] 1/2 , 



(10) 
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which interpolates between the known analytic chain and dimer limits, to the 2x12 lattice 
data. (Except for the points u>'(k = 0) in the secondary band, which showed large finite 
size effects, and which we replaced as argued above by 2u( / k).) For Jj_ — J\\ — J the fitted 
constants were found to be u(0) = 1.890J, uj(it) = 0.507J and c = 1.382J. In Fig. 6 we 
show the triplet dispersion relation which this parametrization predicts for (VO) 2 P 2 7 , 
together with a similar result for the secondary band, using the mean value J± — J\\ — 7.79 
meV and the physical lattice spacing. The primary triplet band extends from 3.9 meV at 
k = 0.813A" 1 to 16 meV at about 0.3A" 1 , and then falls to 15 meV at k — 0. The secondary 
band extends from 7.9 meV at k = to a broad plateau at an energy of about 17-18 meV 
centered on k = 0.813A" 1 . Structure function calculations on the 2x12 lattice suggest that 
the secondary band should appear most clearly near the k = it point (0.813A- 1 ). 

For comparison we quote predictions for the triplet spin-wave dispersion relation in the 
dimer chain model. Of course the lattice spacing a and the direction of the continuous 
momentum variable k are problematical for (VO) 2 P20 7 in the dimer model because there 
is no obvious dimer chain interaction pathway. Since the dimer unit cell has length 2a the 
dispersion relation repeats with period Ak = it /a; this implies that the two different gaps we 
found for the ladder at and it /a are equal in the dimer chain. Another characteristic feature 
of the dimer chain dispersion relation is that it is symmetric about it /2a, due to inversion 
symmetry. For the parameters Ji, J 2 and g found in our susceptibility fits the dimer chain 
model predicts a somewhat larger gap of E gap xs 0.44Jx = 4.9 meV and a bandwidth of 
11 meV. It is interesting that one can apparently distinguish between the dimer chain 
and ladder models by an accurate measurement of the gap alone, using parameters derived 
from susceptibility fits. 

V. Summary and Conclusions 

In this paper we used numerical techniques to study the susceptibility of a Heisenberg an- 
tiferromagnetic spin ladder and a dimerized Heisenberg spin chain. We used exact numerical 
diagonalization to generate all energy eigenvalues and their degeneracies, which were then 
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used to determine x(T) on ladders and dimer chains of up to 16 spins. We presented results 
for a range of temperatures and interaction ratios J±/J\\ (ladder) and J2/J1 (chain). These 
were extrapolated to give bulk limit estimates, which we parametrized using a function with 
five parameters. We fitted the bulk limit x(T) to the susceptibility data for (VO) 2 P20 7 , 
which is a candidate spin ladder system but is known to be accurately described by the dimer 
chain susceptibility. Our best fit to the dimer chain model accurately reproduces previous 
parameter values. Our best fit for the ladder is in slightly better agreement with the data, 
and indicates that (VO) 2 P 2 07 has very similar Jj_ and J\\ values. With these parameters 
we give numerical predictions for the spin-wave excitation gap of (VO) 2 P 2 07 and for other 
properties of the spin-wave dispersion relation. 
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TABLES 

TABLE I. x(T) versus T/Jy for the 2 x 8 ladder model. 



T/J|| 


Oi£ = J±/J\\ 


0.5 


0.7 


0.9 


1.0 


1.1 


0.05 


.001034 


.000434 


.000095 


.000036 


.000012 


0.10 


.024049 


.015862 


.007576 


.004697 


.002730 


0.15 


.052841 


.040907 


.026020 


.019306 


.013692 


0.20 


.071958 


.059399 


.043488 


.035436 


.027946 


0.25 


.086240 


.073189 


.057624 


.049470 


.041489 


0.30 


.098180 


.084851 


.069643 


.061651 


.053662 


0.35 


.108058 


.094958 


.080161 


.072384 


.064535 


0.40 


.115854 


.103404 


.089220 


.081722 


.074105 


0.45 


.121734 


.110128 


.096743 


.089607 


.082308 


0.50 


.125983 


.115240 


.102749 


.096035 


.089124 


0.55 


.128895 


.118955 


.107366 


.101100 


.094615 


0.60 


.130729 


.121508 


.110776 


.104954 


.098906 


0.65 


.131696 


.123113 


.113171 


.107773 


.102151 


0.70 


.131970 


.123957 


.114733 


.109728 


.104508 


0.75 


.131691 


.124192 


.115617 


.110972 


.106127 


0.80 


.130975 


.123942 


.115954 


.111638 


.107136 


0.85 


.129917 


.123308 


.115852 


.111835 


.107648 


0.90 


.128592 


.122373 


.115402 


.111656 


.107754 


0.95 


.127064 


.121205 


.114674 


.111175 


.107534 


1.00 


.125384 


.119856 


.113728 


.110454 


.107051 


1.25 


.115891 


.111691 


.107134 


.104726 


.102238 


1.50 


.106267 


.102995 


.099495 


.097661 


.095775 


1.75 


.097456 


.094849 


.092086 


.090649 


.089174 
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2.00 


.089662 


.087541 


.085311 


.084156 


.082975 


2.25 


.082839 


.081084 


.079250 


.078303 


.077337 


2.50 


.076876 


.075402 


.073867 


.073077 


.072274 


2.75 


.071649 


.070394 


.069092 


.068424 


.067745 


3.00 


.067046 


.065965 


.064848 


.064276 


.063695 


3.25 


.062971 


.062031 


.061062 


.060566 


.060064 


3.50 


.059344 


.058520 


.057671 


.057238 


.056799 


3.75 


.056099 


.055370 


.054621 


.054239 


.053853 


4.00 


.053181 


.052531 


.051866 


.051527 


.051184 
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TABLE II. x(T) versus T/Ji for the N = 16 dimer chain model. 



T/Ji 






ad 


= J2/J1 






0.2 


0.4 


0.6 


0.7 


0.8 


1.0 


0.05 


.000000 


.000001 


.000032 


.000207 


.001380 


.011118 


0.10 


.000616 


.001530 


.006061 


.013706 


.031256 


.073202 


0.15 


.009732 


.015373 


.032088 


.049187 


.073849 


.105906 


0.20 


.035487 


.045476 


.067867 


.084726 


.103142 


.117932 


0.25 


.071654 


.081545 


.100186 


.111536 


.121656 


.124279 


0.30 


.107862 


.114179 


.125146 


.130670 


.134374 


.129372 


0.35 


.138004 


.139715 


.143241 


.144277 


.143674 


.134013 


0.40 


.160397 


.158026 


.155844 


.153834 


.150521 


.138121 


0.45 


.175673 


.170271 


.164219 


.160299 


.155375 


.141511 


0.50 


.185222 


.177821 


.169378 


.164351 


.158550 


.144073 


0.55 


.190459 


.181882 


.172110 


.166518 


.160317 


.145789 


0.60 


.192568 


.183407 


.173028 


.167222 


.160925 


.146707 


0.65 


.192465 


.183123 


.172604 


.166806 


.160603 


.146916 


0.70 


.190827 


.181567 


.171200 


.165544 


.159549 


.146517 


0.75 


.188146 


.179129 


.169092 


.163656 


.157931 


.145614 


0.80 


.184775 


.176099 


.166488 


.161314 


.155889 


.144304 


0.85 


.180966 


.172681 


.163547 


.158652 


.153537 


.142674 


0.90 


.176898 


.169028 


.160386 


.155772 


.150965 


.140797 


0.95 


.172699 


.165246 


.157093 


.152755 


.148245 


.138736 


1.00 


.168456 


.161412 


.153734 


.149660 


.145433 


.136543 


1.25 


.148209 


.142928 


.137255 


.134274 


.131202 


.124795 


1.50 


.130981 


.126962 


.122690 


.120461 


.118172 


.113424 


1.75 


.116819 


.113687 


.110382 


.108667 


.106912 


.103283 


2.00 


.105181 


.102682 


.100061 


.098707 


.097324 


.094474 
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2.25 


.095528 


.093493 


.091370 


.090276 


.089161 


.086869 


2.50 


.087431 


.085744 


.083990 


.083090 


.082173 


.080293 


2.75 


.080559 


.079139 


.077668 


.076914 


.076147 


.074579 


3.00 


.074663 


.073452 


.072201 


.071561 


.070912 


.069584 


3.25 


.069555 


.068511 


.067434 


.066885 


.066327 


.065190 


3.50 


.065090 


.064181 


.063245 


.062768 


.062284 


.061299 


3.75 


.061157 


.060357 


.059537 


.059119 


.058696 


.057834 


4.00 


.057666 


.056959 


.056233 


.055864 


.055491 


.054731 
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